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The divergence theorem of Gauss plays a central role in the derivation of the governing differ- 
ential equations in fluid dynamics, electrodynamics, gravitational fields, and optics. Many of these 
phenomena, particularly shocks and turbulence in fluids, are multi-scale in nature and prone to con- 
tinuous generation of high wave-number modes or small scales. In practical applications, however, 
one is interested in an evolution equation for the large scale quantities without resolving the details 
of the small scales. As a result, there has been a significant effort in developing time-averaged and 
spatially-filtered equations for large scale dynamics from the fully resolved governing differential 
equations. One should realize that by starting from these fully-resolved equations (e.g. the Euler or 
Navier-Stokes equations) to derive an averaged evolution equation one has already taken the limit 
of the wave-numbers approaching infinity with no regards to our observational abilities at such a 
limit. As a result, obtaining the evolution equations for large scale quantities (low wave-numbers) 
by an averaging or filtering process is done after the fact in order to somehow regularize the irregu- 
lar behavior at high wave-numbers. This could explain many of the theoretical and computational 
difficulties with the Euler or Navier-Stokes equations. Here, a rather different approach is proposed. 
The averaging process in implemented before the derivation of the differential form of the transport 
equations. A new observable divergence concept is defined based on fluxes calculated from observ- 
able quantities at a desired averaging scale, a. An observable divergence theorem is then proved 
and applied in the derivation of the observable and regularized transport equations. We further 
show that the application of the observable divergence theorem to incompressible flows results in 
a formal derivation of the inviscid Leray turbulence model first proposed in 1934. Finally, some 
numerical simulations are presented for the observable inviscid compressible Euler equations with 
shocks. It is argued that such a methodology in deriving fluid evolution equations removes many of 
the theoretical and computational difficulties in multi-scale problems such as turbulence and shocks. 
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a. Introduction. For the last several years we have 
been promoting the idea of simultaneous inviscid regular- 
ization of shocks and turbulence in inviscid flows. While 
shocks and turbulence are often treated as separate prob- 
lems p], they share some similar characteristics. The 
primary shared feature between them is a high wavenum- 
ber cascade processes created by the nonlinear convective 
term u • Vu. Traditionally, this high wavenumber irreg- 
ularity is regularized and controlled by the addition of a 
viscous term or a random walk in the Euler equations. 
While successful in describing many complex flow phe- 
nomena, this approach is still facing significant theoret- 
ical and computational challenges in high wavenumbers 
(or high Reynolds number flows) [2 . In this investigation, 
we take a different approach and explore a formal deriva- 
tion of an inviscid regularization of this high wavenumber 
behavior in shocks and turbulence. The essence of the 
technique is in the re-derivation of the divergence theo- 
rem of Gauss for field quantities that are averaged. The 
divergence theorem of Gauss in calculus is revisited in or- 
der to prove an analogous theorem for a control volume 
where the fluxes at the faces of the fluid element are cal- 
culated based on averaged flow quantities over a certain 
length scale as opposed to a fully resolved one. 
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Our work here was motivated by our earlier results in 
inviscid regularization of the ID and multidimensional 
Burgers equations [3J 0] where the nonlinear convective 
term is replaced by u • Vu; u is an averaged velocity. We 
have shown that under certain conditions the solution to 
our regularized Burgers equation exists and is unique at 
all times and it converges to the entropy solution of the 
Burgers equation. Norgard and Mohseni [S] has recently 
extended some of these results to ID compressible Euler 
equations. Supporting numerical simulations were pre- 
sented and it was shown that these equations share the 
same traveling wave solutions and entropy solutions with 
the Euler equations. This manuscript provides a formal 
derivation of our regularized Euler equations from basic 
principles and its extension to multi-dimensional flows. 
It should be noted, while a regular viscous term could be 
included in our analysis without any difficulty that is not 
necessary for our purposes here. 

b. Averaging Kernels. In the following sections we 
are required to calculate the flux of some conserved quan- 
tities which depends on an averaging process. A general 
class of averages for a function / can be defined by the 
convolution operation as / = g * f, where g is the kernel 
of the convolution. In such an averaging there are sev- 
eral guidelines that seem intuitively reasonable; see [3J @]. 
As listed in Table [IJ the averaging kernel should be non- 
negative, decreasing and radially symmetric. Such a filter 
g can be associated with a characteristic wavelength a. 
This parameter is introduced by scaling the filter as such, 
g a = —g (— ), or in Fourier domain <? a (k) = g(ak). Thus 
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TABLE I. Requirements for the averaging kernel. 



Properties 


Mathematical Expression 


Normalized 
Nonnegative 
Decreasing 
Symmetric 


1 9=1 
g(x) > 0, Vx 
|xi| > |x 2 | =>■ g(xi) < gr(x 2 ) 

|xi| = |X 2 => ff(xi) = ff(x 2 ) 



as a becomes smaller, the wavelength where the filter 
exerts influence also become smaller. With this scaling, 
the filter remains normalized, non-negative, decreasing 
and isotropic. Furthermore, as a — ¥ 0, g a approaches 
the Dirac delta function. When convolving g a with /, 
features in / that have length scales less than a will be 
averaged out. 

c. Observability of a Vector Field. Our ability to 
observe and measure (physically or numerically) any fluid 
property is limited by the spatial resolution of our exper- 
imental or numerical technique. For instance, if one uses 
hotwire annemometry for measuring turbulence velocity 
components, the spatial resolution of the measurement is 
always limited by the size of the wire. In fact, what is ac- 
tually measured is an averaged value of the velocity over 
a length scale of the order of the size of the hotwire. Sim- 
ilarly in numerical simulations, one is limited by scales of 
the order of the mesh size. As a result, what we often ob- 
serve is an averaged quantity rather than the mathemat- 
ical limit where the observation volume approaches zero. 
This observation has significant ramifications when one 
calculates the flux of a quantity at the face of a volume 
surrounding a point of interest. Before we investigate 
this, we will revisit the divergence of a vector field and 
the divergence theorem of Gauss. We will then extend 
these concepts to the case of a fluid volume with limited 
observability. 

d. Divergence of a Vector Field. In writing the dif- 
ferential form of the conservation laws for a fluid, one 
considers the flux of the conserved quantities from the 
faces of a fluid volume. These fluxes are then related 
to the divergence of the conserved quantity inside the 
volume. Let Q be a small and closed region centered at 
point xo whose boundary is represented by the orientable 
surface S and its volume is Ay. Let F be a vector func- 
tion in C 1 (r2). We denote the flux of the vector field F 
through the surface S by The divergence of the 
vector field F around the point Xo € is defined as the 
limit of the flux density o/F over £1 around the point Xo, 
that is 
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divF(xo) = lim 



= 1™ v , - 



F-ndS (1) 



S(Q) 



where n is the outer unit normal vector on S. Therefore, 
the divergence of a vector field is basically the volume 
averaged flux of the vector field out of Q; namely flux 
per unit volume. Note that by taking SV — > we have 
already made the assumption that the vector field F is 
observable with infinite resolution. This assumption will 
be challenged later in this manuscript. 

A well-known result from the vector calculus is the 
divergence theorem of Gauss reproduced below [5] . 



Theorem .1 Divergence Theorem of Gauss. Let CI be 
a region with surface boundary S oriented outward. Let 
F(x) G C\Q). Then 



divF dV = 



F • n dS. 



(2) 



e. Divergence of an Averaged Vector Field. In or- 
der to extend the concept of the divergence of a vec- 
tor field, discussed above, to the case of a finite volume, 
where the fluxes are calculated based on averaged quan- 
tities, we will adopt the definition of the divergence as 
defined in equation 0; that is a volume-averaged flux. 
We now follow the standard steps in the proof of the di- 
vergence theorem of Gauss with a more careful examina- 
tion of the flux approximation at the faces of the volume. 
Before this, we first introduce the following definition. 

Definition .1 Observable Divergence. Consider a func- 
tion f and a vector field V. We define the observable 
divergence of a vector field F = /V as 



odiv F = / divV + V • grad/, 



(3) 



where (.) is the averaging process with an observation- 
ability at a length scale a. 

As it will be clear shortly, this quantity represents the 
average flux of the vector field F out of a finite region 
with a representative length scale a defined by the aver- 
aging process (7). For instance, for the Helmholtz filter 
/ = / — a 2 A/ where the representative scale of the finite 
averaging region is a. The observation scale a could be 
imagined to be of similar size to the mesh size of the nu- 
merical simulation or the size of the measurement region 
(e.g., dimension of the wire in hotewire annemometry) in 
experiments. Note that, in a limit where the observable 
scale approaches zero this observable divergence reduces 
to the classical divergence defined in calculus. 

We now prove the main result governing the divergence 
theorem for an observable conservation law. 

Theorem .2 Observable Divergence Theorem. Let be 
a region with surface boundary S oriented outward. Let 
F(x) = /V(x) = (u(x), i>(x), w(x))/ be a vector function 
in C^fi). Then 

odivF=/ / / (/divV + V-grad/) dV = [ jF-ndS 
n J J Jn JJs 

where (7) is defined by the Helmholtz operator. 

Proof. Consider a small closed element of fluid con- 
sisting of a rectangular parallelepiped centered around 
{x,y, z), where the values of / and V are defined, with 
edges parallel to the coordinate axes; see Figure [T] If 
/ and V are continuous the flux at the boundary point 
x s of a finite volume is defined as /(x s ) V(x s ). We now 
start with the divergence definition in equation (JlJ. The 
surface integral ff s F ■ n dS over the rectangular par- 
allelepiped can be regarded as a sum of six terms, one 
for each face of the parallelepiped. Now considering that 
F\ = fu and using the average flux approximation at the 
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FIG. 1. A finite element of fluid at point (x, y, z). 
surface S\ and 5*2, one can write 
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where we used the Taylor expansion and the definition of 
the Helmholtz filter to obtain an approximation for the 
flux on the face of the element [7]. Adding the contribu- 
tions from the rest of the faces for the parallelepiped and 
taking the limit of the volume AV approaching zero, one 
can write 

lim -J— [[ F-ndS 

av^oAV JJs 

= fU X + V,f X + fVy + Vfy + fW Z + Wf Z 

= /V-V + V- V/ = odivP 

F-ndS = /// odivF. ■ 
is JJJn 

One should note that the limit of vanishing volume is 
performed independently and beyond what is observable 
with the averaging scale a. That is the volume of the 
element continuously shrinks until it is not observable 
at the available numerical or experimental resultion. At 
such a limit there is a clear distinction between odivF 
and divF. At the mathematical limit where both the 
volume element AV and the length scale of the averag- 
ing kernel a approch zero independently, the difference 
between the div F and odiv F diminishes. 

/. Observable Conservation Laws. In this section 
we will develop the differential equations that must be 
satisfied by a fluid with an observability limit of a. The 
equations are expected to resolve all relevant dynamical 
quantities up to this observable scale a. As mentioned 
before, this scale is dictated by the resolution of a partic- 
ular numerical simulation or the length scale resolution 
of a particular experimental equipment. For scales below 
a we assume that the medium still acts as a continuum 
and one can take the mathematical limit of AV —> 0. 
However, our observable scale remains at a even as this 
mathematical limit is taken to zero. 

Before the differential form of the observable conserva- 
tion laws for mass, momentum, and energy of a contin- 
uum are presented, we will formulate the general form of 
the observable differential conservation law for a quantity 



per unit volume, /. We assume / satisfies a conservation 
law of the form 

Rate of generation of / = Q (4) 

where Q is the sum of all sources of / inside the volume. 
For instance Q is the total external force exerted on a 
control volume if one considers the conservation of mo- 
mentum inside the volume. The rate of generation of / is 
basically the outflow of / minus the inflow at the bound- 
aries of the volume plus the storage inside the volume. 
The rate of outflow minus the inflow at the boundaries 
was conveniently calculated in the previous section by 
the operator odiv /V where V is the continuum veloc- 
ity. Therefore, the observable conservation of / reduces 
to the following differential equation 



| + 7diw- 



V • grad/ = Q. 



(5) 



/ divV + V • grad/ = 0. 



(6) 



We now show that this equation satisfies the conservation 
of / over a volume. 

Theorem .3 Conservation of Conserved Quantities. 
Consider the observable evolution equation for f 

df_ 
dt 

where f = g* f is the averaged or filtered quantity, and g 
is an averaging kernel satisfying the properties in Table 
[7| It can be shown that the observable evolution equation 
for f conserves the total of f over the whole domain; that 
is q r 

-y/(x)dx = 0. (7) 

Proof. Integrating the equation ^ over the domain and 
using the definition of the averaging, one obtains 



/(y)s(x - y)dy 



V(y) 5 (x-y)dy dx. (8) 



+V/(x) 



Using the identities /V ■ V + V ■ V/ = V ■ (/V) and 
J//(/ vv + v ' V/)tfx = for functions that rapidly 
decay toward the boundary one can easily change the 
order of integration to prove that the right hand side of 
equation ([8| is zero. ■ 
g. Observable Euler Equations. Now, the observ- 
able compressible Euler equations can be easily derived 
by using the conservation of mass (p), momentum (pV), 
and energy {per = p(e + \ V 2 )) to obtain 
dp 



dpV 

~dt 
dpe T 

dt 



dt 

pV divV 



p divV 



- V • gradp 
• grad(pV) 
grad(pe T ) 



0. 



-Vp. 



(9a) 
(9b) 



pex divV + V 
(p divV + V • gradp) + pV ■ S + S, (9c) 



where S is the sum of all external forces acting on the 
control volume or its surface and S is the sum of all en- 
ergy sources inside the control volume. Note that in the 
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limit of a approaching zero, that is when the observable 
scale approaches a mathematical zero, one recovers the 
classical Euler equations. 

In our previous publications [Fj, we presented some 
mathematical and numerical results for the existence and 
uniqueness of the solution in ID Euler along with the 
convergence to the entropy solution when a — > 0. We also 
showed that our ID observable Euler equations support 
the same traveling shock solution as the Euler equations. 

h. Observable Incompressible Euler Equations. For 
the case of constant density, the observable continuity 
equation ( 9a I reduces to the classical case of a divergence 



free velocity held; that is divV = 0. In this case, the 
observable momentum equation ( 9b ) can be simplified to 



<9V — 1 

— + V • gradV = Vp. 

at p 



(10) 



This is exactly the inviscid Leray equation 8J, introduced 
without a formal derivation in 1934. A formal derivation 
of these equations from basic principles is offered here. 

i. Numerical Experiments. Several numerical exper- 
iments with increasing degrees of complexity were per- 
formed in order to numerically evaluate the perfor- 
mance of the observable Euler equations. These include 
benchmark test cases from [9] using the Euler/Navier- 
Stokes equations; ID Sod, ID Shu-Osher, 2D shock- 
vorticity /entropy wave interaction, 3D Taylor-Green vor- 
tex, 3D isotropic compressible turbulence, and 3D shock- 
turbulence interaction. See [10] for the details of our 
numerical experiments. In order to reduce potential nu- 



merical dissipation a pseudo-spectral discretization tech- 
nique is used. The equations are advanced in time with a 
Runge-Kutta-Fehlberg predictor/corrector (RK45). Spa- 
tial derivatives and the inversion of the Helmholtz oper- 
ator were computed in the Fourier domain. The terms 
were converted into the Fourier domain using a Fast 
Fourier Transform, multiplied by the appropriate term 
and then converted back into the physical domain. 

Here we only present the results for the 2D interaction 
of a vorticity/entropy wave with a normal shock located 
at x — 1.57T. Identical initial and boundary conditions 
and similar spatial resolutions as in [5] is used here. Fig- 
ure [2] shows the instantaneous vorticity contours and pro- 
files for two different incident angles. The inflow vortic- 
ity/entropy wave interacts with the shock and changes 
its propagation direction and wavelength. In Figure [2jb) 
instantaneous vorticity profiles for different a are shown 
for the incident angle <f> — 75°. Compared to the solu- 
tions in [3] , the vorticity profiles based on observable Eu- 
ler equations do not show any spurious oscillations across 
the shock. Again note that this is done without any vis- 
cous terms. These results demonstrates the ability of 
the observable Euler equations to simulate 2D flows with 
shocks and resolve the shock without any viscous terms 
or side effects (numerical dissipation) of using one-sided 
schemes. Also note that the thickness of the shock is con- 
trolled by the parameter a, the width of the averaging 
Krenel. For smaller a, higher resolutions are required. 

We expect to employ the observable divergence theo- 
rem in future to other field quantities in order to obtain 
regularized and observable evolution equations in other 
areas beyond the fluid dynamic problems considered here. 
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FIG. 2. Pseudo spectral simulation of the observable Euler 
equation. (Left) Vorticity contours for <j> = 45° at t — 25 
and ki = ki = 1. (Right) Vorticity Profiles at X2 = 0, t=32, 
4> = 75°, and fci = 1 where 4> is the angle of the incoming 
wave and ki is wavelength of the incoming waves. 



